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Abstract. We investigate the global texture model of structure formation in cosmogo- 
nies with non-zero cosmological constant for different values of the Hubble parameter. 
We find that the absence of significant acoustic peaks and little power on large scales 
are robust predictions of these models. However, from a careful comparison with ex- 
periments we conclude that at present we cannot safely reject the model on the grounds 
of published CMB anisotropy data. If bias is close to one on large scales, galaxy corre- 
lation data rules out the models. New, very stringent constraints come from peculiar 
velocities. 

Investigating the large-iV limit, we argue that our main conclusions apply to all 
global O(N) models of structure formation. 



I INTRODUCTION 

Recently, a lot of effort has gone into the determination of cosmological param- 
eters from measurements of cosmic microwave background (CMB) anisotropies, 
especially in view of the two planned satellite experiments MAP and PLANCK [1] . 
However, we believe it is important to be aware of the heavy modeling which en- 
ters these results. In general, simple power law initial spectra for scalar and tensor 
perturbations and vanishing vector perturbations are assumed, as predicted from 
inflation. To reproduce observational data, the composition of the dark matter and 
the cosmological parameters as well as the input spectrum and the scalar to tensor 
ratio are varied [2]. 

In this work, we assume that cosmic structure was induced by scaling seeds. 
We follow the philosophy of a general analysis of scaling seed models motivated in 
Ref. [3]. 

Seeds are an inhomogeneously distributed form of matter (like e.g. topological 
defects) which interacts with the cosmic fluid only gravitationally and which rep- 
resents always a small fraction of the total energy of the universe. They induce 
geometrical perturbations, but their influence on the evolution of the background 



universe can be neglected. Furthermore, in first order perturbation theory, seeds 
evolve according to the unperturbed spacetime geometry. 

Here, we mainly investigate the models of structure formation with global tex- 
ture. This models (for VL matter = 1) show discrepancies with the observed intermedi- 
ate scale CMB anisotropies and with the galaxy power spectrum on large scales [4] . 
Recently it has been argued that the addition of a cosmological constant leads to 
better agreement with data for the cosmic string model of structure formation [5]. 
We analyze this question for the texture model, by using ab initio simulation of 
cosmic texture as described in Ref. [6]. We determine the CMB anisotropies, the 
dark matter power spectrum and the bulk velocities for these models. We also 
compare our results with the large- N limit of global O(N) models, and we discuss 
briefly which type of parameter changes in the 2-point functions of the seeds may 
lead to better agreement with data. 

We find that the absence of significant acoustic peaks in the CMB anisotropy 
spectrum is a robust result for global texture as well as for the large- N limit for 
all choices of cosmological parameters investigated. Furthermore, the dark matter 
power spectrum on large scales, A ~ 20/i~ 1 Mpc, is substantially lower than the 
measured galaxy power spectrum. 

However, comparing our CMB anisotropy spectra with present data, we cannot 
safely reject the model. On large angular scales, the CMB spectrum is in quite good 
agreement with the COBE data set, while on smaller scales we find a significant 
disagreement only with the Saskatoon experiment. Furthermore, for non-satellite 
experiments foreground contamination remains a serious problem due to the limited 
sky and frequency coverage. 

The dark matter power spectra are clearly too low on large scales, but in view of 
the unresolved biasing problem, we feel reluctant to rule out the models on these 
grounds. A much clearer rejection may come from the bulk velocity on large scales. 
Our prediction is by a factor 3 to 5 lower than the POTENT result on large scales. 

Since global texture and the large- N limit lead to very similar results, we conclude 
that all global O(N) models of structure formation for the cosmogonies investigated 
in this work are ruled out if the bulk velocity on scales of 50/i _1 Mpc is around 
300km/s or if the CMB primordial anisotropies power spectrum really shows a 
structure of peaks on sub-degree angular scales. 

The formalism used for our calculations is not presented here but can be found 
in [7]. There, we also explain in detail the eigenvector expansion which allows to 
calculate the CMB anisotropies and matter power spectra in models with seeds 
from the two point functions of the seeds alone. Section 2 is devoted to a brief 
description of the numerical simulations. In Section 3 we analyze our results and 
in Section 4 we draw some conclusions. 

Notation: We always work in a spatially flat Friedmann universe. The metric is 
given by 

ds 2 = a(t) 2 (dt 2 - 5ijdx l dx j ) , 
where t denotes conformal time. 



Greek indices denote spacetime coordinates (0 to 3) whereas Latin ones run from 
1 to 3. Three dimensional vectors are denoted by bold face characters. 



II THE NUMERICAL SIMULATIONS 



Like in [6], we consider a spontaneously broken iV-component scalar field with 
O(N) symmetry. We use the a-model approximation, i.e., the equation of motion 



where (5 is the rescaled field (5 = <p/r]. 

We do not solve the equation of motion directly, but use a discretized version of 
the action [8]: 



where A is a Lagrange multiplier which fixes the field to the vacuum manifold (this 
corresponds to an infinite Higgs mass). Tests have shown that this formalism agrees 
well with the complementary approach of using the equation of motion of a scalar 
field with Mexican hat potential and setting the inverse mass of the particle to the 
smallest scale that can be resolved in the simulation (typically of the order of 10 -35 
GeV), but tends to give better energy momentum conservation. 

As we cannot trace the field evolution from the unbroken phase through the 
phase transition due to the limited dynamical range, we choose initially a random 
field at a comoving time t = 2Ax. Different grid points are uncorrected at all 
earlier times [9]. 

To determine the CMB anisotropies and dark matter power spectra, we need 
to calculate the unequal time correlators (UTC) of the energy momentum tensor 
or the gravitational field induced by the seeds. In Figs 1, 2, 3, we present as an 
example one of the scalar UTC's and the scalar time correlators, r = t'/t = 1, as 
well as the values at z = (k 2 tt r ) = as a function of r. Since there are two scalar 
degrees of freedom, the two Bardeen potentials $ and \1/ there are three correlators: 



The UTCs are obtained numerically as functions of the variables k, t and t c with 
t > t c and t c fixed. They are then linearly interpolated to the required range. 
We construct a hermitian 100 x 100 matrix in kt and kt', with the values of kt 
chosen on a linear scale to maximize the information content, < kt < x max . The 



np-(p-np)p = o, 



(i) 



S = J d 4 xa 2 (t) l -d^ ■ 9>/3 + ± ((3 2 - l) 



(2) 




(3) 
(4) 
(5) 




FIGURE 1. The two point correlation function C u (z,r) = fc 4 v / tt 7 ($ s (k, t)$*(k, t')) is shown. 
Panel (a) represents the result from numerical simulations of the texture model; panel (b) shows 
the large- N limit. For fixed r the correlator is constant for z < 1 and then decays. Note also the 
symmetry under r — » 1/r. 
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FIGURE 2. The correlators Cy(z, 1) are shown. The solid, dashed and dotted lines represent 
C22 i Cn and \C\2\ respectively. Panel (a) is obtained from numerical simulations of the texture 
model and panel (b) shows the large- TV limit. A striking difference is that the large- iV value for 
|C 12 1 is relatively well approximated by the perfectly coherent result \J\CnC22\ while the texture 
curve for | C12 1 lies nearly a factor 10 lower. 




choice of a linear scale ensures good convergence of the sum of the eigenvectors 
after diagonalization (see [7]), but still retains enough data points in the critical 
region, 0(x) = 1, where the correlators start to decay. In practice we choose as 
the endpoint x max of the range sampled by the simulation the value at which the 
correlator decays by about two orders of magnitude, typically x max ~ 40. The 
eigenvectors that are fed into a Boltzmann code are then interpolated using cubic 
splines with the condition v n {kt) — ► for kt 3> a; max . 

We use several methods to test the accuracy of the simulation: energy momentum 
conservation of the defects code is found to be better than 10% on all scales larger 
than about 4 grid units, as is seen in Fig. 4. A comparison with the exact spherically 
symmetric solution in non-expanding space [8] shows very good agreement. 

The resulting CMB spectrum on Sachs Wolfe scales is consistent with the line of 
sight integration of Ref. [6]. Furthermore, the overall shape and amplitude of the 
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FIGURE 3. The correlators Cy(0, r) are shown in the same line styles as in Fig. 2, but for z = 
as function of r — t' /t. The stronger decoherence of the texture model is even more evident here. 
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FIGURE 4. Energy momentum conservation of our numerical simulations is shown. The lines 
represent the sum of the terms which has to vanish if energy (solid) respectively momentum 
(dashed) is conserved, divided by the sum of the absolute value of these terms. The abscissa 
indicates the wavelength of the perturbation as fraction of the size of the entire grid. 



unequal time correlators are quite similar to those found in the analytic large- N 
approximation [10,11,3] (see Figs. 1,2 and 3). The main difference of the large- 
N approximation is that there the field evolution, Eq. (1), is approximated by a 
linear equation. The non-linearities in the large- N seeds which are due solely to the 
energy momentum tensor being quadratic in the fields, are much weaker than in the 
texture model where the field evolution itself is non-linear. Therefore, decoherence 
which is a purely non-linear effect, is expected to be much weaker in the large- N 
limit. This is actually the main difference between the two models as can be seen 
in Fig. 5. 
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FIGURE 5. Top panel: the f w = model. Bottom panel: The Ce power spectrum is shown 
for the large- N limit (bold line) and for the texture model. The main difference is clearly that 
the large- N curve shows some acoustic oscillations which are nearly entirely washed out in the 
texture case. 



Ill RESULTS AND COMPARISON WITH DATA 



A CMB anisotropies 



The Cgs for the 'standard' global texture model are shown in Fig. 6 (bottom 
panel), where a comparison of the full result with the totally coherent approxima- 
tion (see [7]) is presented. 




FIGURE 6. The Ce power spectrum for the texture scenario is shown in the perfectly coherent 
approximation (top panel) and in the full eigenfunction expansion. Even in the coherent approx- 
imation, the acoustic peaks are not higher than the Sachs Wolfe plateau. Decoherence (see [7]) 
just washes out the structure but does not significantly damp the peaks. 
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TABLE 1. The value of the normalization constant e and the fluctuation ampli- 
tude erg are given for the different models considered. The error in e comes from a 
best fit normalization to the full CMB data set. Cosmological parameters which 
are not indicated are identical in all models or given by Oo = £l c .dm + + ^fc = 1 ■ 
We consider only spatially flat models with Of, = 0.05 and a helium fraction of 
23%. The parameter choice indicated in the top line is referred to as standard 
texture model in the text. 

Vector and tensor modes are found to be of the same order as the scalar com- 
ponent at COBE-scales. For the 'standard' texture model we obtain : Cyp : 
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FIGURE 7. The Ci power spectrum is shown for different values of cosmological parameters. 
In the top panel we choose SI a = 0, Qcdm — 0.95, fif, = 0.05 and vary h. In the bottom panel 
we fix h = 0.5, fib = 0.05 and vary We only consider spatially flat universes, f^o = 1. 

Cfjp ~ 0.9 : 1.0 : 0.3, in good agreement with the predictions of Refs. [4,12,21] and 
[3]. Due to tensor and vector contributions, even assuming perfect coherence (see 
Fig. 6, top panel), the total power spectrum does not increase from large to small 
scales. Decoherence leads to smoothing of oscillations in the power spectrum at 
small scales and the final power spectrum has a smooth shape with a broad, low 
isocurvature 'hump' at i ~ 100 and a small residual of the first acoustic peak at 
I ~ 350. There is no structure of peaks at small scales. The power spectrum is 
well fitted by the following fourth-order polynomial in x = logi: 



The effect of decoherence is less important for the large- N model, where oscil- 
lations and peaks are still visible (see Fig 5, bottom panel). This is due to the 
fact that the non-linearity of the large- N limit is only in the quadratic energy mo- 
mentum tensor. The scalar field evolution is linear in this limit [10], in contrast to 
the N = 4 texture model. Since decoherence is inherently due to non-linearities, 
we expect it to be stronger for lower values of N. COBE normalization leads to 
e = (0.92 ± 0.1) 10" 5 for the large- N limit. 

In Fig. 7 we plot the global texture Ce power spectrum for different choices of 
cosmological parameters. The variation of parameters leads to similar effects like 
in the inflationary case, but with smaller amplitude. At small scales (£ > 200), the 
Cgs decrease with increasing H and they increase when a cosmological constant 
Q\ = 1 — Q m is introduced. Nonetheless, the amplitude of the anisotropy power 
spectrum at high is remains in all cases on the same level like the one at low is, 
without showing the peak found in inflationary models. The absence of acoustic 
peaks is a stable prediction of global O(N) models. The models are normalized to 
the full CMB data set, which leads to slightly larger values of the normalization 
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1.5 - 2.6a; + 3.3a; 2 - 1.4a; 3 + 0.17a; 4 . 



(6) 



parameter e = AttGt] 2 than pure COBE normalization. In Table 1 we give the 
cosmological parameters and the value of e for the models shown in Fig. 7. 

In order to compare our results with current experimental data, we have selected a 
set of 31 different anisotropy detections obtained by different experiments, or by the 
same experiment with different window functions and/or at different frequencies. 
Theoretical predictions and data of CMB anisotropies are usually compared by 
plotting the theoretical Ce curve along with the CMB measurements converted to 
band power estimates. We do this in the top panel of Fig. 8. The data points raise 
from large to smaller scales, in contrast to the theoretical predictions of the model. 
This fashion of presenting the data is surely correct, but lacks informations about 
the uncertainties in the theoretical model. Therefore we also compare the detected 
mean square anisotropy, A( Exp > and the experimental \-a error, Y,( Exp > 2 , directly 
with the corresponding theoretical mean square anisotropy, given by 

A (Tk )= 1 Y,(2e + l)C e We , (7) 

47T £ 

where the window function Wg contains all experimental details (chop, modulation, 
beam, etc.) of the experiment. 

The theoretical error in principle depends on the statistics of the perturbations. 
If the distribution is Gaussian, one can associate a sample/cosmic variance 

where / represents the fraction of the sky sampled by a given experiment. 

Deviation from Gaussianity leads to an enhancement of this variance, which can 
be as large as a factor of 7 (see [14]). Even if the perturbations are close to Gaussian 
(which has been found by simulations on large scales [6,15]), the Q's, which are 
the squares of Gaussian variables, are non-Gaussian. This effect is, however only 
relevant for relatively low is. Keeping this caveat in mind, and missing a more 
precise alternative, we nevertheless indicate the minimal, Gaussian error calculated 
according to (8). We add a 30% error from the CMB normalization. The numerical 
seeds are assumed to be about 10% accurate. 

In Table 2, the detected mean square anisotropy, A^ Exp \ with the experimental 
1-er error are listed for each experiment of our data set. The corresponding sky 
coverage is also indicated. In Fig. 8 we plot these data points, together with the 
theoretical predictions for a texture model with h = 0.5 and Q\ = 0. 

We find that, apart from the COBE quadrupole, only the Saskatoon experiment 
disagrees significantly, more than la, with our model. But also this disagreement 
is below 3a and thus not sufficient to rule out the model. In the last column of 
Table 2 we indicate 

X 2 = (A (T7l) - A (Exp) ) 2 /(£ {T/l)2 + ^ Exp)2 ) 



(8) 
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FIGURE 8. The Cp, spectrum obtained in the standard texture model is compared with data. 
In the top panel experimental results and the theoretical curve are shown as functions of I. In the 
two lower panels we indicate the value of each of the 31 experimental data points with l-cr error 
bars and the corresponding theoretical value with its uncertainty. The experiments corresponding 
to a given number are given in Table 2. In the middle panel the 8 COBE data points are shown. 
In the bottom panel other experiments are presented. 
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TABLE 2. The CMB anisotropy detections used in our analysis. The 3., 4. and 5. column denote the 
value of the anisotropy and the upper and lower 1-cr errors respectively. The references are: Tegmark and 
Hamilton 1997 [16]; de Bernardis et al. 1994 [17]; Cheng et al. 1994 [18]; Cheng et al. 1996 [19]; Cheng et 
al. 1997 [20]; Tanaka et al. 1996 [22]; Gutierrez et al. 1997 [23]; Gunderscn et al. 1993 [24]; Dragovan et 
al. 1993 [25]; Masi et al 1996 [26]; Netterfield et al. 1996 [27]; Scott et al. 1997 [28]. 



for the j-th experiment, where the theoretical model is the standard texture model 
with Qa = and h = 0.5. The major discrepancy between data and theory 
comes from the COBE quadrupole. Leaving away the quadrupole, which can be 
contaminated and leads to a similar x 2 also for inflationary models, the data agrees 
quite well with the model, with the exception of three Saskatoon data points. 
Making a rough chi-square analysis, we obtain (excluding the quadrupole) a value 
X 2 = J2j Xj ~ 30 for a total of 30 data points and one constraint. An absolutely 
reasonable value, but one should take into account that the experimental data 
points which we are considering are not fully independent. The regions of sky 
sampled by the Saskatoon and MSAM or COBE and Tenerife, for instance, overlap. 
Nonetheless, even reducing the degrees of freedom of our analysis to iV = 25, our 
X 2 is still in the range (N — 1) ± y/2(N — 1) ~ 24 ± 7 and hence still compatible 
with the data. 

This shows that even assuming Gaussian statistics, the models are not convinc- 
ingly ruled out from present CMB data. There is however one caveat in this analy- 
sis: A chi-square test is not sensitive to the sign of the discrepancy between theory 
and experiment. For our models the theoretical curve is systematically lower than 
the experiments. For example, whenever the discrepancy between theory and data 
is larger than 0.5a, which happens with nearly half of the data points (13), in all 
cases except for the COBE quadrupole, the theoretical value is smaller than the 
data. If smaller and larger are equally likely, the probability to have 12 or more 
equal signs is 2(13 + l)/2 13 ~ 3.4 x 10~ 3 . This indicates that either the model is 
too low or that the data points are systematically too high. The number 0.003 can 
however not be taken seriously, because we can easily change it by increasing our 
normalization on a moderate cost of x 2 ■ 

B Matter distribution 

In Table 1 we show the expected variance of the total mass fluctuation an in a 
ball of radius R = 8/i~ 1 Mpc, for different choices of cosmological parameters. We 
find <t 8 = (0.44 ± 0.07)h (the error coming from the CMB normalization) for a 
flat model without cosmological constant, in agreement with the results of Ref. [4]. 
From the observed cluster abundance, one infers a$ = (0.50 ± 0.04)f2~ a5 [29] and 
a$ = 0.59tojg [30]. These results, which are obtained with the Press- Schechter 
formula, assume Gaussian statistics. We thus have to take them with a grain of 
salt, since we do not know how non-Gaussian fluctuations on cluster scales are in 
the texture model. According to Ref. [31], the Hubble constant lies in the interval 
h ~ 0.73 ± 0.06 ± 0.08. Hence, in a flat CDM cosmology, taking into account 
the uncertainty of the Hubble constant, the texture scenario predicts a reasonably 
consistent value of cr 8 . 

As already noticed in Refs. [21] and [4], unbiased global texture models are 
unable to reproduce the power of galaxy clustering at very large scales, ~ 20h~ 1 
Mpc. In order to quantify this discrepancy we compare our prediction of the linear 



matter power spectrum with the results from a number of infrared ( [32], [33]) 
and optically-selected ( [34], [35]) galaxy redshift surveys, and with the real-space 
power spectrum inferred from the APM photometric sample ( [36]) (see Fig. 10). 
Here, cosmological parameters have important effects on the shape and amplitude 
of the matter power spectrum. Increasing the Hubble constant shifts the peak of 
the power spectrum to smaller scales (in units of h/Mpc), while the inclusion of a 
cosmological constant enhances large scale power. 

We consider a set of models in - h space, with linear bias [37] as additional 
parameter. In Table 3 we report for each survey and for each model the best value 
of the bias parameter obtained by ^-minimization. We also indicate the value of 
X 2 (not divided by the number of data points). The data points and the theoretical 
predictions are plotted in Fig. 21. Our bias parameter strongly depends on the 
data considered. This is not surprising, since also the catalogs are biased relative 
to each other. 

Models without cosmological constant and with h ~ 0.8 only require a relatively 
modest bias b ~ 1.3 — 3. But for these models the shape of the power spectrum is 
wrong as can be seen from the value of x 2 which is much too large. The bias factor 
is in agreement with our prediction for For example, our best fit for the IRAS 
data, for h ~ 0.8 is b ~ 1.3. With al RAS = (0.69 ±0.05), this gives a 8 ~ 0.48 ±0.04, 
compatible with the direct computation 

Whether IRAS galaxies are biased is still under debate. Published values for 
the /3 parameter, defined as (3 = Q°' 6 /b, for IRAS galaxies, range between (3j = 
0.9lg; 2 5 [38] and /3j = 0.5 ±0.1 [39]. Biasing of IRAS galaxies is also suggested 
by measurements of bias in the optical band. For example, Ref. [40] finds f3 = 
0.40 ± 0.12, in marginal agreement with [41], which obtains (3 Q = 0.35 ± 0.1. A 
bias for IRAS galaxies is not only possible but even preferred in flat global texture 
models. 

But also with bias, our models are in significant contradiction with the shape 
of the power spectrum at large scales. As the values of x 2 m Table 3 and Fig. 10 
clearly indicate, the models are inconsistent with the shape of the IRAS power 
spectrum, and they can be rejected with a high confidence level. The APM data 
which has the smallest error bars is the most stringent evidence against texture 
models. Nonetheless, these data points are not measured in redshift space but they 
come from a de-projection of a 2 — D catalog into 3-D space. This might introduce 
systematic errors and thus the errors of APM may be underestimated. 

Models with a cosmological constant agree much better with the shape of the 
observed power spectra, the value of x 2 being low for all except the APM data. But 
the values of the bias factors are extremely high for these models. For example, 
IRAS galaxies should have a bias b ~ 3 — 6, resulting in as < 0.25, and in a f3j < 0.2 
which is too small, even allowing for big variances due to non-Gaussian statistics. 

The power spectra for the large- N limit and for the coherent approximation are 
typically a factor 2 to 3 higher (see Fig. 10), and the biasing problem is alleviated 
for these cases. For Ha = we find a 8 = 0.57h for the large- N limit and a$ = 0.94/i 
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16 


19 


CfA2-SSRS2 130 Mpc 


0.5 


0.4 


5.6 


5 


19 


CfA2-SSRS2 130 Mpc 


0.5 


0.8 


11.1 


4 


19 


LCRS 


0.5 


0.0 


3.0 


71 


19 


LCRS 


0.8 


0.0 


1.8 


96 


19 


LCRS 


1.0 


0.0 


1.6 


108 


19 


LCRS 


0.5 


0.4 


3.7 


33 


19 


LCRS 


0.5 


0.8 


8.7 


40 


19 


IRAS 


0.5 


0.0 


2.3 


102 


11 


IRAS 


0.8 


0.0 


1.3 


131 


11 


IRAS 


1.0 


0.0 


1.3 


140 


11 


IRAS 


0.5 


0.4 


2.8 


70 


11 


IRAS 


0.5 


0.8 


6.3 


9 


11 


IRAS 1.2 Jy 


0.5 


0.0 


4.2 


56 


29 


IRAS 1.2 Jy 


0.8 


0.0 


2.9 


92 


29 


IRAS 1.2 Jy 


1.0 


0.0 


2.9 


99 


29 


IRAS 1.2 Jy 


0.5 


0.4 


4.3 


39 


29 


IRAS 1.2 Jy 


0.5 


0.8 


6.7 


28 


29 


APM 


0.5 


0.0 


3.3 


1350 


29 


APM 


0.8 


0.0 


1.8 


1500 


29 


APM 


1.0 


0.0 


1.7 


1466 


29 


APM 


0.5 


0.4 


3.5 


1461 


29 


APM 


0.5 


0.8 


6.2 


1500 


29 


QDOT 


0.5 


0.0 


4.3 


32 


19 


QDOT 


0.8 


0.0 


2.9 


44 


19 


QDOT 


1.0 


0.0 


2.9 


46 


19 


QDOT 


0.5 


0.4 


4.3 


25 


19 


QDOT 


0.5 


0.8 


7.3 


14 


19 



TABLE 3. Analysis of the matter power spectrum. In the first column the catalog 
is indicated. Cols. 2 and 3 specify the model parameters. In cols. 4 and 5 we give 
the bias parameter inferred by \ 2 minimization as well as the value of x 2 . Col. 6 
shows the number of 'independent' data points assumed in the analysis. 
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FIGURE 9. The dark matter power spectrum for the texture model (solid line) is compared 
with the coherent approximation (short dashed) and the large- N limit (long dashed). The spectra 
are COBE normalized and the cosmological parameters are f?A = , h = 0.5. 

for the coherent approximation. This is no surprise since only one source function, 
the analog of the Newtonian potential, seeds dark matter fluctuations and 
thus the coherence always enhances the unequal time correlator. The dark matter 
Greens function is not oscillating, so this enhancement translates directly into the 
power spectrum. 

Models which are anti-coherent in the sense defined in Section IID reduce power 
on Sachs- Wolfe scales and enhance the power in the dark matter. Anti-coherent 
scaling seeds are thus the most promising candidates which may cure some of the 
problems of global O(N) models. 

The simple analysis carried out here does not take into account the effects of non- 
linearities and redshift distortions. Redshift distortions in the texture case should 
be less important than in the inflationary case since the peculiar velocities are rather 
low (see next paragraph). Non-linearities typically set in at k > 0.5/iMpc -1 and 
should not have a big effect on our main conclusions which come from much larger 
scales. Inclusion of these corrections will result in more small-scale power and in 
a broadening of the spectra, which even enhances the conflict between models and 
data. Furthermore, variations of other cosmological parameters, like the addition 
of massive neutrinos, hot dark matter, which is not considered here, will result in 
a change of the spectrum on small scales but will not resolve the discrepancy at 
large scales. 

Nonetheless, scale dependent biasing may exist and lead to a non-trivial relation 
between the calculated dark matter power spectrum and the observed galaxy power 
spectrum. We are thus very reluctant to rule out the model by comparing two in 
principle different things, the relation of which is far from understood. Therefore 
we would prefer to reject the models on the basis of peculiar velocity data, which 
is more difficult to measure but most certainly not biased. 
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FIGURE 10. Matter Power spectrum: comparison between data and theory. References are in 
the text. Data set courtesy of M. S. Vogeley [44]. 
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<j v (R) 


A„ 


h = 0.5 


h = 1.0 


Ov = 0.8 


10 


494 


170 


145 


205 


86 


20 


475 


160 


100 


134 


78 


30 


413 


150 


80 


98 


70 


40 


369 


150 


67 


78 


65 


50 


325 


140 


57 


65 


61 


60 


300 


140 


50 


56 


57 



TABLE 4. Bulk velocities: Observational data from [42] and theoretical predictions. A„ esti- 
mates the observational uncertainty. The uncertainties on the theoretical predictions are around 
<~ 30%. The models Ha — with h = 0.5 and h = 1 as well as CIa = 0.8, h = 0.5 are investigated. 



C Bulk velocities 

To get a better handle on the missing power on 20 to 100/i _1 Mpc, we investigate 
the velocity power spectrum which is not plagued by biasing problems. The as- 
sumption that galaxies are fair tracers of the velocity field seems to us much better 
justified, than to assume that they are fair tracers of the mass density. We therefore 
test our models against peculiar velocity data. We use the data by Ref. [42] which 
gives the bulk flow 

°1{R) = / P(k)W(kR)dk , (9) 

in spheres of radii R — 10 to 60/i _1 Mpc. These data are derived after recon- 
structing the 3— dimensional velocity field with the POTENT method (see [42] and 
references therein). 

As we can see from Table 4, the COBE normalized texture model predicts too 
low velocities on large scales when compared with POTENT results. Recent mea- 
surements of the bulk flow lead to somewhat lower estimates like <J V (R) ~ (230 ±90) 
at R — 60/i _1 Mpc ( [43]), but still a discrepancy of about a factor of 2 in the best 
case remains. 

Including a cosmological constant helps at large scales, but decreases the veloci- 
ties on small scales. 

If the observational bulk velocity data is indeed reliable (there are some doubts 
about this [45]), all global O(N) models are ruled out. 



IV CONCLUSIONS 

We have determined CMB anisotropies and other power spectra of linear per- 
turbations in models with global O(N) symmetry which contain global monopoles 
and texture. Our main results can be summarized as follows: 

• Global O(N) models predict a flat spectrum (Harrison-Zeldovich) of CMB 
anisotropies on large scales which is in good agreement with the COBE results. 



Models with vanishing cosmological constant and a large value of the Hubble 
parameter give cr 8 ~ 0.4 to 0.5 which is reasonable. 

• Independent of cosmological parameters, these models do not exhibit pro- 
nounced acoustic peaks in the CMB power spectrum. 

• The dark matter power spectrum from global O(N) models with f2 A = has 
reasonable amplitude but does not agree in its shape with the galaxy power 
spectrum, especially on very large scales > 20/i _1 Mpc. 

• Models with considerable cosmological constant agree relatively well with the 
shape of the galaxy power spectrum, but need very high bias b ~ 4 — 6 even 
with respect to IRAS galaxies. 

• The large scale bulk velocities are by a factor of about 3 to 5 smaller than the 
values inferred from [42]. 

In view of the still considerable errors in the CMB data (see Fig. 8), and the biasing 
problem for the dark matter power spectrum, we consider the last argument as the 
most convincing one to rule out global O(N) models. Even if velocity data is quite 
uncertain, observations agree that bulk velocities on the scale of 50/i _1 Mpc are 
substantially larger than the (50 - 70)km/s obtained in texture models. 

However, all our constraints have been obtained assuming Gaussian statistics. 
We know that global defect models are non-Gaussian, but we have not investigated 
how severely this influences the above conclusions. Such a study, which we plan 
for the future, requires detailed maps of fluctuations, the resolution of which is 
always limited by computational resources. Generically we can just say that non- 
Gaussianity can only weaken the above constraints. 

Our results naturally lead to the question whether all scaling seed models are 
ruled out by present data. The main problem of the O(N) model is the missing 
power at intermediate scales, I ~ 300 - 500 or R ~ (20 - 100)/i~ 1 Mpc. We 
have briefly investigated whether this problem can be mitigated in a scaling seed 
model without vector and tensor perturbations. In this case, also scalar anisotropic 
stresses are reduced by causality requirements (see Ref. [3]), and perturbations on 
superhorizon scaled are compensated (see [46]). For simplicity, we analyze a model 
with purely scalar perturbations and no anisotropic stresses at all, f n = 0. The 
seed function $ s is taken from the texture model (numerical simulations) and we set 
\J/ S = —& s . The resulting CMB anisotropy spectrum is shown in Fig. 5, top panel. 
A smeared out acoustic peak with an amplitude of about 2.2 does indeed appear 
in this model. This is mainly due to fluctuations on large scales being smaller, as 
is also evident from the higher value of e = (2.2 ± 0.2) x 10~ 5 . But also here, the 
dark matter density fluctuations and bulk velocities are substantially lower than 
observed galaxy density fluctuations or the POTENT bulk flows. 

Clearly, this simple example is not sufficient and a more thorough analysis of 
generic scaling seed models is presently under investigation. So far it is just clear 
that contributions from vector and tensor perturbations are severely restricted. 
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